Tsunami waveform inversion by numerical finite-elements Green?s functions
نویسندگان
چکیده
During the last few years, the steady increase in the quantity and quality of the data concerning tsunamis has led to an increasing interest in the inversion problem for tsunami data. This work addresses the usually ill-posed problem of the hydrodynamical inversion of tsunami tide-gage records to infer the initial sea perturbation. We use an inversion method for which the data space consists of a given number of waveforms and the model parameter space is represented by the values of the initial water elevation field at a given number of points. The forward model, i.e. the calculation of the synthetic tide-gage records from an initial water elevation field, is based on the linear shallow water equations and is simply solved by applying the appropriate Green’s functions to the known initial state. The inversion of tidegage records to determine the initial state results in the least square inversion of a rectangular system of linear equations. When the inversions are unconstrained, we found that in order to attain good results, the dimension of the data space has to be much larger than that of the model space parameter. We also show that a large number of waveforms is not sufficient to ensure a good inversion if the corresponding stations do not have a good azimuthal coverage with respect to source directivity. To improve the inversions we use the available a priori information on the source, generally coming from the inversion of seismological data. In this paper we show how to implement very common information about a tsunamigenic seismic source, i.e. the earthquake source region, as a set of spatial constraints. The results are very satisfactory, since even a rough localisation of the source enables us to invert correctly the initial elevation field.
منابع مشابه
Uncertainties in the 2004 Sumatra–Andaman source through nonlinear stochastic inversion of tsunami waves
Numerical inversions for earthquake source parameters from tsunami wave data usually incorporate subjective elements to stabilize the search. In addition, noisy and possibly insufficient data result in instability and non-uniqueness in most deterministic inversions, which are barely acknowledged. Here, we employ the satellite altimetry data for the 2004 Sumatra-Andaman tsunami event to invert t...
متن کاملRapid Detection of Large Earthquakes Using Quasi-Finite-Source Green’s Functions in Moment Tensor Analysis
Rapid procedures (i.e. within 5 to 15 minutes) for earthquake and tsunami early warnings focus on earthquake location, depth, magnitude, and slowness, but initial tsunami early warnings are more often issued without knowing the mechanism of the earthquake. Kawakatsu (1998) proposed using a limited number of stations to automatically detect, locate and determine the source parameters of earthqua...
متن کاملDiscretized Adjoint State Time and Frequency Domain Full Waveform Inversion: A Comparative Study
This study derives the discretized adjoint states full waveform inversion (FWI) in both time and frequency domains based on the Lagrange multiplier method. To achieve this, we applied adjoint state inversion on the discretized wave equation in both time domain and frequency domain. Besides, in this article, we introduce reliability tests to show that the inversion is performing as it should be ...
متن کاملEffects of dispersion in tsunami Green's functions and implications for joint inversion with seismic and geodetic data: a case study of the 2010 Mentawai MW 7.8 earthquake
Tsunami observations play an important role in resolving offshore earthquake slip distributions. Non-dispersive models are often used with an initial static sea-surface pulse derived from seafloor deformation in computation of tsunami Green's functions. We compare this conventional approach with more advanced techniques, which use Green's functions computed by a dispersive model with an initial...
متن کاملSensitivity of the adjoint method in the inversion of tsunami source parameters
This paper tests a methodology for tsunami waveform inversion, based on the adjoint method. The method is designed to perform the direct optimization of the tsunami fault parameters, from tide-gauge data, imposing strong geophysical constrains to the inverted solutions, leading to a substantial enhancement of the signal-to-noise ratio, when compared with the classical technique based on Green’s...
متن کامل